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Abstract 

A geometric definition for a magnetic charge of Abelian monopoles in SU(N) lat- 
tice gauge theories with Higgs fields is presented. The corresponding local monopole 
number defined for almost all field configurations does not require gauge fixing and is 
stable against small perturbations. Its topological content is that of a 3-cochain. A 
detailed prescription for calculating the local monopole number is worked out. 

Our method generalizes a magnetic charge definition previously invented by Phillips 
and Stone for SU(2). 
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1 Introduction 



The investigation of the complicated vacuum structure of quantized Yang-Mills theories re- 
quires the use of non-perturbative methods. The lattice regularization is the most convenient 
in this respect. Among other applications it allows to identify those excitations (monopoles, 
instantons, etc.) which are expected to dominate the Euclidean path integral at large scales 
and to explain important features like chiral symmetry breaking and quark confinement in 
QCD. 

Monopoles in which we take interest in this paper form the basis for the dual super- 
conductor scenario of confinement P], |2|. According to this scenario the condensation of 
(anti-) monopoles causes the color-electric flux to be squeezed into string-like flux tubes 
(dual Meissner effect). The consequence is a linearly rising potential between static quarks. 

Whether extended non-Abelian 't Hooft-Polyakov-like monopoles (see @) or singular 
Abelian Dirac-like monopoles play the major role in the confinement phenomenon has not 
been settled yet. The majority of investigations concentrates on the latter kind of excitations. 
These can be identified by Abelian projection after an appropriate gauge has been fixed [||]. 
Up till now the maximally Abelian gauge was considered to be the most favourable gauge 
HI H 01 (f° r a recent review see e.g. f|). For this gauge the resulting U{1) degrees of 
freedom as well as the corresponding Abelian monopoles were shown to reproduce almost 
completely the string tension of the full theory (so-called Abelian and monopole dominance, 



respectively) P, ITO 



Other gauges - defined for instance by diagonalizing gauge observables as parallel trans- 
porters along closed paths (e.g. the Polyakov loop) which transform like Higgs fields in the 
adjoint representation - have been studied, too |pH . Based on appropriately defined order 
parameters according to a construction of Frohlich and Marchetti [|I^] the condensation of 
Abelian monopoles in the confinement phase was demonstrated to occur widely independent 
of the gauge chosen [|13[ . 

But a general gauge independent proof for the validity of the dual superconductor picture 
is still missing. 

Therefore, it would be worthwhile to have a gauge-independent and - at the same time - 
geometric lattice definition of the monopole charge. This, in general, requires the existence 
of a Higgs field in the adjoint representation. Magnetic fluxes and the magnetic charge, re- 
spectively, can easily be obtained by projecting the non-Abelian plaquette variables onto the 
local colour direction of the Higgs field. However, the naive discretization does not provide 
any geometric meaning of the resulting charge. Nevertheless, corresponding investigations 
of magnetic fluxes have provided useful insight into the phase structure of the 4d SU{2) 
Georgi-Glashow model |14| . 



A first geometrical definition for the monopole number has been given in |J| [16| . The 
lattice gauge field together with the Higgs field determines local maps from the boundary 
of cubes into the coset space SU(2)/U(1). The winding numbers of these maps provide 
the local integer-valued monopole charge, for each cube, or alternatively speaking, a closed 
magnetic current on the dual lattice. In our present work we are going to generalize the 
definition to a SU (N) gauge-Higgs theory. The Higgs fields employed are considered to be 
CP iV_1 -valued fields. We shall present all the details necessary for developing a computer 
code for the calculation of the winding numbers from the local lattice field data. 
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The resulting magnetic charge definition can be employed for explorations of the phase 
structure of SU (N) gauge-Higgs models. But, it might be worthwhile also for the application 
to pure SU(N) gauge theories, when an auxiliary Higgs field is determined as proposed for 
the recently invented Laplacian gauge fixing fT7|j . 

Recently, in the context of investigations of the electro-weak phase transition for the 
SU(2) theory another version of an integer-valued, gauge independent lattice magnetic 



charge definition was invented |{L8 |. It remains an interesting task to work out the inter- 



relation between the different definitions, if there is any. 

Section 2 will provide the general description and the necessary notations. In Section 3 
we present the computation of the monopole number. Some detailed proofs are shifted into 
the Appendix. In Section 4 we shall draw the conclusions. 



2 Description of the model and notation 

In this work we want to consider a model consisting of an SU(N) lattice gauge field U = 
{Uij} coupled to a Higgs field $ = {</>;} transforming with respect to the fundamental 
representation. i,j, . . . denote the lattice sites and pairs (ij) the corresponding links of a 
4-dimensional simplicial lattice. Without loss of generality for our construction of the local 
monopole number we work with a Higgs field of fixed length, \\fa\\ = 1. In the following we 
denote by Uijk... = UijUjk ■ ■ ■ parallel transporters longer than one link. 

As the phase of the vector fa e will be irrelevant for the monopole charge, it is 
sufficient to identify the corresponding equivalence class The latter can be represented 
by the N x N projection matrix $j = fa (g) <p* projecting on the subspace spanned by the 
vector fa. In this sense the Higgs field takes values in a CP^" 1 manifold. 

The action of a complex invertible matrix V on elements of C~P N ~ l will be denoted by 
V, e.g. V * $ means the projector Vcf) £g> (V</>)*/||\^0|| 2 or the class [Vfa. For the lattice 
we assume an ordering of the lattice sites and denote the different kinds of simplices (sites, 
links, triangles, tetrahedra, 4-tetrahedra) by ordered tuples of vertices such as (ij . . . k) . 
If we consider a site i of a given simplex (ij . . . k) as the coordinate origin, a cube c^-...^ 
parametrized by coordinates < Sj, . . . , < 1 can be defined. 

3 Definition of the local monopole number 

Our aim in this section is to define for a given tetrahedron A = (ijkl) a map from the 
boundary of the 3-dimensional cube c/f^n to CP^ -1 . This map is constructed from the link 
and Higgs variables at the particular tetrahedron A. 

Since ^(CP^ -1 ) = Z, the homotopy class of such a map is represented by an integer 
riA- This integer will be referred to as the local monopole number at the particular simplex. 
It will be shown to have the following properties: 

(a) Gauge invariance: ua is independent of a particular choice of gauge. 

(b) Local conservation: 
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The assignment A i— > n& represents a cochain in the sense that 

(Sn) T = (-)a^a = 0, 

AedT 

the sign depending on the relative orientation of A in the boundary of a 4-tetrahedron 
T. If /i is the link dual to A we may equivalently say that M M = n& is a conserved 
(divergence free) magnetic current on the dual lattice. 

(c) Geometric significance: 

The collection of the Xa's define a section of a SU(N)/U(N — l)-bundle over the 
2-skeleton of the lattice. This section may be regarded as an interpolation of the Higgs 
field which is as horizontal as possible with respect to the lattice gauge field. 

The maps Xa are worked out stepwise for higher simplices. In order to state their 
relevant properties, we need GL(X)-valued parallel transport functions Vuj_ k i){ s j, . . . , s k ) 
as introduced in [fEBI. These functions essentially measure to what extent the lattice gauge 
field fails to be flat. They are an important ingredient in the construction of SU(N) bundles 
associated to lattice gauge fields. In our case, the relevant formulae read 

%> = u a> 

V(ijk) = SjU ijk + (1 — Sj)U ik , 

V(ijki) = SjSkUijM + (1 - Sj)s k Uiki + Sj(l - s k )Uiji + (1 - Sj)(l - s k )Uu, 

where < Sj, s k < 1. Let us assume for simplicity that the simplex A in question is (01 ... k). 
Then Xa is assumed to meet the following requirements: 

X(0...k)(- ■ ■ > S r — 1, • • •) = V( ...r)(si, . . .) * X( r ... fc )(s r , . . .), (1) 

X(o...k)(- ■ ■ , s r = 0, . . .) = X{p...f...k)(si,...,$ r ,...), (2) 

and likewise for all other simplices, with the obvious replacements. It follows from the results 
in Jl6| that such a collection of maps defines a section over the 2-skeleton in a bundle with 
principal fiber CP^ 1 . It is also seen that the local winding number of this section at a 
tetrahedron A = (ijkl) is just the winding number of X^i) : dc^ki) — > CP^ -1 . Let us now 
explicitly construct this map. 

Requirements (HD and (fj) suggest to define the X&s inductively for simplices of increasing 
dimension. This is done in four steps: 

Step 0: As the Xa's essentially interpolate the Higgs field, we set X^ = $j for the 
lattice sites. 

Step 1: On the far end Sj = 1 of the interval cuj\, we must have, by requirement ([!]) 

X {ij)( S 3 = !) = %> * x {j)- 

On the base Sj = 0, X^j) must be X^, by requirement (H). We thus interpolate geodesically 
in CP 7 ^ 1 for intermediate values of Sj. 

Step 2: On the far sides Sj = 1 or Sk = 1 of the square c/^-m, we must have, by 
requirement (P 



(A) X( ijk )(sj = 1, s k ) = Vft) * X(jk)(sk), 

(B) X {ijk ) (sj, s k = 1) = V(ijk) (.Sj) * X {k ) . 
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Again, on the base Sj = s k = 0, X^ k ) must be X@, by requirement (Q). For intermediate 
values we interpolate geodesically from the the far sides to the base. 

Step 3: In the last step we proceed as before. Again, on the far faces Sj = 1, s k = 1 or 
si = 1 of the cube c^jki), we must have, by requirement ([I]) 

(-0 X {ijkl)( s j = 1) Sfe, Si) = Vft) * X( jkl )(s k , Si), 

(II) X(i jkl )(sj, s k = 1, si) = V{ijk)(sj) * X( kl )(si), 
(III) X(ij k i)(sj, s k , si = 1) = V(ij k i)(sj, s k ) * X(i). 

As above, on the base Sj = s k = si = 0, X^jki) = Xq, by requirement @. Unlike in the 
previous cases, we cannot interpolate from the the far sides to the base since there might be 
an obstruction. Nevertheless, we might do so on the boundary of our cube. This completes 
our construction. Requirements ([Tp and (Q) are more or less obviously satisfied. 

In the next chapter we shall see how the winding number of Xuj k i\ may be calculated. 



4 Calculation of the monopole number 

We have to determine the class of X&. Take a simplex A = (0123). The boundary of the 
unit cube ca with base (0) is the union of six faces C k ,C k (k = 1, 2, 3), where 

Cfc = {(si, s 2 , s 3 ) : s k = 0}, Cl = {(s u s 2 , s 3 ) : s k = 1}. 



The monopole number may be computed as an intersection number |T9|j . 

n A = X A (dc A ) : E 1 , 

E denoting the (2N — 4)-dimensional space of all projectors perpendicular to some rank 
one projector E. The intersection number does not depend on the particular projector E 
chosen. To facilitate the calculations, we take it to be $o- The strategy is to compute 
the intersection number separately for the six different faces of the surface of our cube. As 
geodesic interpolation is needed for the definition of X A , let us collect some facts about 
the geometry of complex projective spaces (see also |2(| [2T, 0]). Two points P,Q in the 



complex projective space may be joined by a unique shortest geodesic if and only if PQ ^ 0. 
Furthermore, for such points it holds: 

Lemma 4.1 The unique shortest geodesic between two points P, Q intersects E 1 - if and only 
if Retr[PQE] < and lmtr[PQE] = 0. Moreover, such an arc can intersect E only once. 

We shall also need the following fact about triangles. 

Lemma 4.2 Let P,Q,R be three points such that there is a unique shortest geodesic between 
any two of them. Let A be the surface spanned by all geodesies between P and the geodesic 
arc QR. This surface intersects E 1 - if and only if the equation lmip(t) = has a solution 
< t < 1 satisfying Reip(t) < 0. Here 

(p(t) = t 2 tr[QR] 1/2 tr[ERP}+ (3) 
(1 - t) 2 tr[QR] l 'Hr[EQP] +t(l-t) tr[E(QR + RQ)P}. 
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The sign of such an intersection is given by 



sign Im 



' dtp® ' 
dt 



Lemma |4.1| and [4.2| are proven in the appendix. Let us first turn our attention to the faces 
Cj[?. With our particular choice of E the following lemma holds. 

Lemma 4.3 Except for a measure zero set of configurations (U, $) , the surfaces X( i23) (C£) 
do not intersect with $q . 

The proof of this lemma is given in the appendix. It follows that we only have to consider 
the faces C\ y C\,C\, schematicly shown in Fig. [|as quadrilaterals (I), (II), (III) respectively. 
On these faces, the map X/qi^) has the form given in Eqs. (I), (II), (HI). The monopole 
number will be the sum of the intersection numbers obtained in each case. 



P^U 01 **, 



Q=U 012 *cD 2 /L\ 


J\ 

IBj\s=U 013 *<I>3 




R=u 0123 ** 3 \ 


/ 11 

m- « 


111 \ 

» * 



X=U Q2 *<£ 2 



W=U *<J> 

VV-U„ 23 M'j 



T=U n 



Figure 1: The image of the surfaces C\ as treated in cases (I), (II), (HI). 



Case (I): The image of C\ is described by Eq. (I). As illustrated in Fig. |2|, it is 
naturally divided into two triangles A,B, according to the structure of Eqs. (A) and (B). 
Let us first consider triangle A. The corners of this triangle are the points P = Uqi * 3>i, Q = 
Uqiz * $2, R = ^0123 * $3- The surface enclosed by these three points is swept out by the set 
of all geodesies between P and the geodesic arc QR, its orientation being fixed by the order 
P,Q,R of the vertices, see Fig. [|. This is the situation described in Lemma [4.2| , except for 
degenerate cases and we use this result to calculate the intersection number of A with $q\ 
Inserting the above expressions for P, Q, R into Eq. (§), we have to consider 



<flA(t) = ^(0o|^O12303)(^12303|0l)(^O10l|0o)|(02|^2303)| 



(4) 
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+ (1 -t) 2 (0 O |f/oi202)(f/l20 2 |0l)(f/oi0l|0o)K02|f/ 2 303)| 
+ t(l -t)((0 O |^O1202)(0 2 |^2303)(^12303|0l)(^O10l|0o) 
+ (0o|^O12303)(^"2303|0 2 )(^ r i202|0l)(^O10l|0o)), 

where the standard notation for the scalar product in has been used. 



f/013*3>3 




{/012*4>2 



Figure 2: The image of C\ according to Eq. (I). 



In order to find the intersections, we must look for solutions < t < 1 of Im (pj^ (t) = 
for which Re<piA{t) < 0. According to Lemma |4~2| such solutions will contribute with the 



sign 



sign Im 



d<piA(t) 
dt 



to the intersection number. 

Triangle B is swept out by geodesies from the dotted line Fig. ||to the point P = t/oi*$i. 
This arc is parametrized by VVn V/i23)(s2) * $3. Hence by Lemma we have to find points 
< s 2 < 1 such that Im(^/ B (s 2 ) = and Re ^75(52) < 0, where ip IB is the function 

tr [x (0) (V (0 i) * X {1) )(V {m (s 2 ) * X {3) ] 
Explicitly, in terms of the parallel transporters and Higgs fields, this may be written as 

<£Ib{s 2 ) = sl(<f)o\Uoi4>l) (4>l\Ui234>3) (Uoi23<p3\4>o) (5) 

+ (1 - 5 2 ) 2 (0 O |^O10l)(0l|f/l303)(f/oi303|0o) 

+ S 2 (l - S 2 )(0o|^O10l)((0l|^12303)(^O1303|0o) 

+ (0l|^1303)(^O12303|0o)), 

which is explicitly seen to be gauge invariant (we have suppressed a positive overall factor 
which obviously does not affect our considerations). Values of s 2 with the above properties 
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may thus be found by elementary means. In the appendix, case a) we show that they make 
the following contribution to the intersection number, i.e. have the following sign for the 
relative orientation: 



sign Im 



Case (II): The image of C\ is a square. It is obtained by 'moving down' the geodesic 
arc PQ to the arc XW, see Fig. |j. From the parametrization of this surface, Eq. (II) one 
concludes that it will intersect the manifold if and only if 



(V< 12)(si) *X(23>(s 3 ))X(o> = 0, 
(V(210>(si) * X (0) )X ( 23)(S 3 ) = 0, 



(6) 



for some < Si, S3 < 1. From Eq. (^) we see that this will happen precisely if the geodesic arc 
parametrised by X^tyiss) intersects the manifold (V{2io)(si)*$o) ± for some value < s\ < 1. 
The arc has endpoints X< 2 ) and V( 23 ) * -^(3)- We may apply Lemma [O] to conclude that 
for such an intersection we must have lmip n (si) = and Reipn( s i) < 0, where tpu is the 
function 



tr 



^(2>(V(23> * X (;i ))(V(2W){si) * X(o)) 

Expressing everything in terms of parallel transporters and Higgs fields, we find 

tPll(Sl) = si (02 1 ^2303) (^012303 1 <t>0) (00 1 U 012 (j) 2 ) 



(7) 



,2\2 



(02 1 ^2303) (^02303 1 0o) (00 |t^O202> 



+ (1 

+ - Sl)(0 2 |?/23 03)((^O123 03|0o)(0o|^O202) 

+ (^O23 03|0o)(0o|^"oi202)), 

and this is directly seen to be gauge invariant. Again as in case (I), we have omitted an 
irrelevant positive overall factor. It is shown in the appendix, case ft) that such an intersection 
contributes with the sign 



sign Im 



dsi 



Case (III): From Eq. (Ill) we can see that the image of C3 is a square, parametrized 
by V(oi23) (si, s 2 ) * X(3), where < Si,s 2 < 1. This square will intersect the manifold $q 
precisely for values si, s 2 such that 



X{ ){V( Q123 )(si,S 2 ) *X( 3 )] 



0. 



We have to find such values and determine their contribution to the intersection number. 
Writing everything in terms of parallel transporters and Higgs fields, the above condition 
reads (pni(si,s 2 ) = 0, where 



</?j//(si, s 2 ) 



+ 



^o|^O12303)SlS 2 + (00 1 ^01303) Si (1 - S 2 ) 
^o|^O2303)(l - Sx)s 2 + (0o|f^O303)(l ~ Si) (1 



(8) 



S2 
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and this expression is manifestly gauge invariant. Condition (|8]) constitutes a system of two 
real equations in two real variables, which can easily be solved. In the appendix, case 7) it 
is demonstrated that any point (si,S2) satisfying Eq. (§) contributes the value 



sign Im 



d(p ni (si,s 2 ) difni(si,s 2 ) 



dsi 



ds<2 



to the intersection number. 



5 Conclusions and outlook 

In this work we have given a geometric definition of local monopole numbers for an SU(N) 
lattice gauge-Higgs system. The definition yields a closed magnetic current which is stable 
under perturbations of the configuration in question. It is explicitly worked out how to 
calculate the monopole number for each tetrahedron of a simplicial 4-dimensional lattice. 
On the practical side, one has to solve four quadratic equations in one variable for each 
such tetrahedron. The coefficients of these equations are given by traces of parallel trans- 
porters and Higgs variables. The prescription is expicitly gauge invariant and can be directly 
implemented in a computer code. 

While the construction gives an integer local monopole number for almost all field con- 
figurations, the question is on what kind of monopole the construction actually triggers. In 
the SU(2) case for the Georgi-Glashow model it has been shown |15| that the construction 



is able to detect a discretized 't Hooft-Polyakov monopole, provided the configuration is 
smooth enough. This question remains to be discussed for N > 2. On the other hand, for 
rough lattice fields as they are generated in a Monte-Carlo simulation, one has to expect 
nonvanishing monopole charges which do not correspond to physical excitations but are lat- 
tice artefacts (dislocations). It is then necessary to find out whether these configurations 
become sufficiently suppressed in the continuum limit. 



Appendix 

Points in CP^ -1 are identified with equivalence classes [a] = Ca (lower case letters) or 
projectors A (capital letters) projecting on the subspace spanned by a. Two points [p], [q] 
are joined by the unique shortest geodesic [(1 — t)(p\q)p + t\{p\q)\ 2 q], provided {p\q) 7^ (see 
12T 



Proof of Lemma 



From the form of the geodesies it is clear that the geodesic arc PQ will intersect E iff 
(1 — t)(e\p)(p\q) + t\(p\q)\ 2 (e\q) = for some < t < 1. This is easily seen to be equivalent 
to 

Re[(e\p)(p\q)(q\e)] < 0, lm[(e\p) (p\q) (q\e)] = 0, 

proving the lemma. 
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Proof of Lemma JJi : 

We parametrize the geodesic arc QR by [p'(t)\ as above. Now the triangle A intersects 
E iff some arc P'(t)P does. According to lemma |4.1| this can happen if and only if 
Im tr [PP'(t)E] = and Re tr [PP'(t)E] < 0. Inserting the expression for P'(t), we see 
that these are just the conditions given in the lemma. The statement for the sign of the 
intersection is obtained in case a. A) below. 



Proof of Lemma 

The image under the map X^ i23) of the faces C° is obtained by coning the boundary of Fig. 
|] to the point E := $o- Let A be any point in this boundary. In order that the geodesic arc 
AE intersects E L , it is necessary that Retr AE < 0, by Lemma P . But this is easily seen 
to imply AE = 0, i.e. A e E x . So the image of the faces C° intersects E L if and only if the 
boundary of the triangle in Fig. |2| does. On dimensional grounds this can happen only for 
a measure zero set of configurations. 

We now prove the various claims about intersection numbers made above. Let us briefly recall 
how the relative orientation of two oriented, embedded transversally intersecting manifolds 
is defined. Suppose X,Y G Z are oriented embedded manifolds, intersecting transversally 
at a point p. The union of two oriented frames in T P X and T P Y then gives a frame in T P Z 
whose orientation may be compared to the orientation of Z. The sign X : Y is defined 
to be ±1 according to whether these orientations coincide resp. do not coincide. Going 
over to the case at hand, we identify with the real vector space H 2N and Euclidean 
inner product (-, ■) = Re( ■ , ■ ). The tangent space T^CP^ -1 is canonically identified 
with {x G = H 2N \ Ax = 0}. Each such tangent space may be equipped with the (real) 
(2N — 2)-form *(aAia), defining thus a canonical orientation of the complex projective space. 
Here, V denotes the operation of taking the Hodge dual of an alternating form with help of 
the above Euclidean inner product. In the same fashion, the tangent space TaE ± , EA = is 
identified with {x e C N = R 2N \ Ax = Ex = 0} and the (real) (2iV-4)-form *(e AieAaAia) 
defines an orientation. Suppose we are given a surface in projective space parametrized by 
[cr(si,s 2 )]. With the above identifications, the tangent space at a point [er ] of this surface 
is given by the real linear span of x%, X2, where 



Xi 



da_ 

dsi 



da_ 



and the 2-form x\ Axi defines an orientation. Now let [o"o] be also in E- 1 , i.e. an intersection 
point. Then by definition, (o~o\xi) = (o~o\e) = and one calculates 



Xi A X2 A *(e A ie A do A iao) = Im 
Hence the sign of this intersection point is given by 

sign Im 



da \ lda_ 
dsi \ds 2 



* (<r A io- ) 



da\ da_ 

dsi / \ ds 2 



(9) 



We are now ready to prove the claims made above about the intersection numbers. The 
maps given in Eqs. (/), (//), (HI) define a surface in CP^ -1 which is schematically drawn 
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in Fig. 0. The arrows define an orientation of this surface and the intersection numbers are 
defined with respect to this orientation and the canonical orientations given above. 
Case a. A): By definition the surface I A is parametrized by [c(s, £)], where 



a(t, s) = (1 - s)((l>(t)\UoM(f>(t) + s|(0(t)|f/ O i0i)| 2 f/oi0i, 

<j){t) = (l-t)(02|f/2303)f/oi202 + t|(02|^230 3 )| 2 f/oi2303- 



The orientation of this surface is given by the order of its vertices P, Q, R, defined above. It 
is easily seen that this means that this orientation coincides with 77 A £, where 



Now using Eq. @ and the definition of <piA in Eq. @, we see that the sign of an intersection 
is given by 



Case a.B): The surface IB is a triangle, parametrized by [cr(s 2 , S3)], where 

v(s 2 ,S 3 ) = (! - S 3)(^(123>(s 2 )03|0l) ^OlV{i23)(-S 2 )03 + -S3KV(i23)(s2)03|0l)| 2 ^0101- 

With the above notations, £2 A £3 is the orientation of this surface. Setting e = <po, us- 
ing Eq. (|9]) and the definition of <piB in Eq. (|5|), it is straightforward to calculate that an 
intersection has sign 



Case (3): Surface II is a square, parametrized by [c(si, S3)], the orientation of this piece 
of surface being given by X\ A x 3 . 

Cr(Si,S 3 ) = (1 - S3)(02|^2303)V(oi2)(Sl)02 + S 3 |(02|^2303)| 2 V(pj3) 0l) ^2303- 

Now setting e = (fi , using Eq. (^) and the definition of <pu in Eq. (|7|), we get the result after 
a short computation. 

Case 7): Surface III is a square, parametrized by [cr(s 1 ,s 2 )] ) where X\ A x 2 is its 
orientation and 




sign Im 



d<PiA{s 2 ) 
ds 2 



sign Im 



dipi B (s 2 ) 
ds 2 



0-(Si,S 2 ) = V(O123>(S1,S2)03- 



Now setting e = (f>o, it follows immediately from the definition of (pm in Eq. (|8|) and from 
Eq. (H) that the sign of the intersection is given by 



sign Im 



dtfnijsu s 2 ) dyinjsx, s 2 ) 
dsi ds 2 
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